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The upper estimate of the percolation threshold of the Bernoulli random field on the hexagonal 
lattice is found. It is done on the basis of the cluster decomposition. Each term of the decomposition 
is estimated using the number estimate of cycles on the hexagonal lattice which represent external 
borders of possible finite clusters containing the fixed lattice vertex. 
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^ ! 1. Introduction. Generally, the percolation theory studies random sub- 

CN . sets of the infinite set where the filtering relative to the inclusion and the 
supplementary connectedness relation has been defined jT]. In particular, such 



Ph. 



(N 



Ph 



^ ; a theory may be arisen in noncompact topological space where filtering is de- 
fined by sequences of compact spaces when their union coincides with the total 
space [2] . The problem of existence of the random realization with connected 
noncompact components represents the main interest. In the case when the 
probability of such an event is positive, they say that the percolation is present 
c5 ! in the random set. The calculation of the percolation probability represents 
the difficult mathematical problem even in most simple mathematical struc- 
^ tures pointed out. Therefore, they resort usually to computer experiments 

> ; for the problem solving when the percolation theory is applied and there are 
some considerable achievements in this direction (see, for example, [3]). In this 
^ ' work, we study the pointed out problem from the mathematical point of view. 
0> ■ Usually, due to the extreme complexity of the problem, such investigations are 
^ ; connected with the study of elementary noncompact spaces such as integer 
^ ' lattices Z^, (i = 2, 3 with the definite connectedness relation. The definition of 
this relation transforms integer lattices into those mathematical objects which 
are called periodic graphs [4j. Percolation on periodic graphs is the subject of 
the discrete percolation theory. However, even for periodic graphs, the main 
problem of percolation theory are resisted to mathematical processing only for 
random sets generated by the Bernoulli field {c(x);x G Z'^} when the proba- 
bility distribution defined by the unique parameter c = Pr{c(x) = 1}. In this 
work, it is found the upper estimate of the so-called percolation threshold c* 
in the case of the two-dimensional uniform periodic graph which is called the 
hexagonal lattice. The percolation probability Q{c) is differed from zero at 
c > c^. Our estimate is done by the well-known approach which we name the 
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cluster decomposition (see also [5], [B]). To estimate the value c^, we find the 
upper estimate of the number of finite clusters on the hexagonal lattice which 
contain the fixed lattice vertex. 

2. The percolation theory problem on the hexagonal lattice. First 
of all, we introduce some geometrical objects modeling crystal lattices. After 
that we will set the problem of the discrete percolation theory on such math- 
ematical structures. In connection with the purpose of the present work, we 
consider only two-dimensional lattices. 

We name the infinite set V in the periodic one if there is the pair (ei, e2) 
of not coUinear vectors in (the parallelogram of periods) such that the 
relation V = V+niei+n2e2 takes place for any G Z, G {1, 2}. We name the 
periodic set in the crystal lattice if it consists of isolated points. The crystal 
lattice admits the disjunctive decomposition V = I J {Vq + niei + n2e2} 



where the finite set Vq is called the crystal cell. If the number of points in Vq is 
minimal among all crystal cells admissible in V, then such a cell is called the 
elementary one. 

Since only the topological structure of the set is important if the percolation 
of random field is studied then, for the formulation of percolation theory prob- 
lem, it is convenient to use the concept of the periodic graph and its immersion 
M into defining the connectedness on the crystal lattice. 

Definition! [4]. Let A = (V, be the infinite nondirectional loop- 
free graph where V is the set of vertexes and $ is the set of some two-element 
subsets {x, y} C V (edges of the graph). This graph is called the periodic 
one with the dimensionality 2 if it admits such an immersion M into M? when 
the image MV is the crystal lattice in and the image M$ of the set $ is 
invariant relative to translations with the parallelogram of periods (ei, e2), 



(ni,n2) G such that the set $o = {{x,y} G $ : My G Fq, x G is finite. 

If {x, y} G X, y G y, then such vertexes are called the adjacent ones 
and we designate the adjacency relation between them by means of x0y. 

Further, we does not distinguish vertexes of the graph and their images 
obtained by the immersion M. We does not distinguish also immersions in M? 
of the graph A differing from each other. Thus, we consider that the vertex 
set V of the graph coincides with I J {Z^ + x} and the property (1) of its 
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M<I) + niei + n2e2 = M$ 
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periodicity is written down in the form $ = $ + niei + 71262. The adjacency 
relation of the periodic graph {V, $) is completely defined by the set $0 since 
the $ admits the disjunctive decomposition $ = |^ {$0 + ^i^i + ^262}. 

(ni,n2)GZ2 

In this connection, we name the set $0 the adjacency one. In terms of 
crystal physics, it defines some "nearest neighbors" on the crystal lattice of 
vertexes being contained in the fixed elementary crystal cell Vq. 

The infinite periodic two-dimensional graph A is called the hexagonal lat- 
tice if its elementary cell Vq contains two vertexes. Besides, at the choice of 
the period parallelogram (2ei,ei + 3e2/2) defined by basis vectors ei,e2 in 
R^, it is possible to put Vq = {xi = ei + 62/2, X2 = 2ei + 62} and $0 = 
{(xi,xi -62); (xi,xi + 62/2 - ei); (xi,X2); (x2,X2 + e2); (x2,X2 + ei -62/2)}. 
The hexagonal lattice {V,^) is shown on the left-hand side of Fig.l where 
the periodic structure have been represented by the dotted line on the right- 
hand side of it. It is formed by shifts of basis vectors of the plane using the 
decomposition on elementary cells which are parallelograms imposed on the 
lattice. 




Figure 1: The hexagonal lattice. 

Let us introduced into consideration the Bernoulli random field {c(x) ; x G 
V} with the concentration c = Pr{c(x) = 1} on the graph {V,^). Here- 
inafter, the tilde which is put over any mathematical object designates that 
it is random. Each random realization c(x), x G V of the field defines the 
set W = {k : c(x) = 1} that we name the configuration. Then, the total set 
{c(x) ; X G y } of random realizations together with the probability distribution 
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on them defines the random set on V. Its probabihty distribution is induced 
by the probabihty distribution of the field {c(x);x £ V}. Namely, for each 
finite subset M C V oi lattice vertexes, the probability of their filling in the 
random configuration W is defined by the formula Pr{M C W} = c'^L 

The adjacency relation ^ induces the connectedness for each random config- 
uration W using the concept of the way on the graph (V, The sequence of 
vertexes {xq, Xi..., x^) chosen from the configuration W is called the connected 
way a with the length n if Xi^x^+i, z = 0,l...,n— 1. The way is called the sim- 
ple one if x^ 7^ Xj in the sequence pointed out at all values i 7^ j, i,j = 1, n 
of indexes and, accordingly, it is called the cycle if the coincidence of vertexes 
in the sequence (xq, Xi..., x^) takes place only at i = 0, j = n. We name the 
vertex pair {x, y} the connected one on W if {x, y} C W and there is a simple 
way a = (x, Xi..., x„_i, y) on this configuration. The connectedness of all ver- 
tex pairs is the equivalence relation. Therefore, any random configuration W 

is broken up the family 9Jl[l^] = {Wj;j G N} of not intersected and connected 

00 

sets, W = which are called clusters. Each cluster consists of vertexes 

i=i 

connected among themselves and any two vertexes being taken from different 
clusters are not connected. We designate the cluster of the family 50t[VK] which 
contains the vertex x G ^ by means of VF(x). If the vertex x is not contained 
in the configuration W we consider that VF(x) = 0. 

The following random function a(x) describes the percolation property of 
the random field {c(z), z &V}, 

n(^) = i ^' if 1^(^)1 = 

1 0; if |iy(x)| <oo 

where the designation | • | = Card(-) is hereinafter used. In the first case, there 
is an infinite simple way a(x) = {xf^i G N+) where x = xq and x^ G VF(x), 
i G N+. In the second case, such a way is absent. On the basis of the random 
function a(x), we may define the probability of the percolation from the fixed 
vertex y G for the field {c(z); z £ V}. They say that the percolation from 
this vertex takes place if the probability Q{c) = Pr{a(y) = 1} is positive. This 
probability does not depend on the vertex y if the periodic graph is uniform 
(for each pair x, y G there are such two immersions Mi, M2 in that their 
images coincide with each other and Mix = M2y). In the present work, we are 
interesting of the value = inf{c : Q{c) > 0} which is called the percolation 
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threshold. 

3. Finite clusters on the hexagonal lattice. Following [Ij, we introduce 
the concept of the external border of the finite cluster W{x.). For this aim, we 
build another periodic graph A* = {V, $*) on the set V which is called the 
conjugate one to the graph A. The adjacency relation on the graph A* is 
introduced as it is shown on Fig. 2 where all vertexes being 0*-adjacent with 
the fixed vertex are pointed out. They are numbered clockwise. This mean 
that vertexes 

{x2+e2,X2+ei+3e2/2,X2+e2 + 2ei,X2 + 2ei,X2+ei-e2/2,X2+ei-3e2/2, 

X2-2e2,X2-ei-3e2/2,X2-ei-e2/2,X2-2ei,X2-2ei + e2,X2-ei + 3e2/2} 

are 0*-adjacent with the vertex X2 on Fig.l at the used immersion of the 
hexagonal lattice. 




Figure 2: The adjacency relation of the vertex 0. 

The adjacency relation 0* as well as the relation generates the connected- 
ness for each random configuration W. It breaks up both the configuration W 
and the configuration V \ being additional to it into some connected sets. 

Definition2 [4J. Let W be the finite cluster. The set dW is called 
the external border of W on the configuration WifWcW and dW consists 
of all vertexes z E Z'^ \ W having the following properties. 

1. For each z, there is the vertex y E W such that z0y. 

2. For the pointed out vertex z, there exists the infinite (p-way a on the 
graph (y, a CiW = which is begun from this vertex and it is unique in 
the intersection a fl {x ^ : x0y, y G W}. 

The following statement is valid. 
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Theorem 1. Let W{'k) be the finite cluster containing the vertex x. G V. 
Then VK(x) has the nonempty finite external border dW{'x.) which possesses 
the following properties. 

1. The set dW{x.) is the simple cycle on the graph (V, 

2. The cycle dW{'K) surrounds the vertex x. at the periodic immersion {V, $*) 
into R2. 

Although this statement is obvious, nevertheless, its full proof uses the so- 
called Jordan topological theorem [1] . 

The way simplicity property is inconvenient to express analytically. There- 
fore, the estimate proposed below in the item 6 is based on the revealing of 
suitable sufficient conditions which are guaranteed that the fixed way forms 
the cluster external border. Using these conditions, we will find the number 
estimate of external borders of all possible finite clusters containing the given 
vertex. 

4. The cluster decomposition on 7?. Let A be the family of finite 
clusters W containing the vertex of the hexagonal lattice. For each such a 
nonempty cluster W, we define the random event 74(iy) = {W : G W, W G 
M[W] , W{0) = W}. The probability of this event is equal to 

FT{A{W)} = c\^\l-cf^\. (2) 

According to Theorem 1, for each cluster W of the family A, there is the sim- 
ple 0*-cycle beloning to the family B = {7 = dW; W G A}. For each 0*-cycle 
7 G B, we define the event 5(7) = {M : G M, W{0) G m[W],dW{0) = 7} 
which is represented in the form of the following finite union of disjoint events 

Bh)= U A{W). (3) 

WeA:dW^j 

The probability ^(7) = Ft{B(^)} is equal to 

p(7)= MAm}= Yl ci^'(i-c)'"^' 

according to (2) and (3). 

It is valid the disjoint decomposition {a(0) = 0} = A{W). The family 

WeA 

A is decomposed on some nonintersecting classes consisting of clusters joined 
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by the following property. We refer clusters W E A to the same class if they 
have the coinciding external border. Therefore, it is valid the transformation 

u...=u ( u 

WeA 7eB \WGA:dW=j J 

Further, on the basis of (3), it follows that {a(0) = 0} = | |J ^(7) U {0 ^ 

\7€B / 

W}. Thus, noticing that Pr{0 ^ ly} = 1 - c and Pr{a(0) = 0} = 1 - (5(c), 
we come to the following statement. 

Theorem 2. The probability 1 — Q{c) is represented by the cluster 
decomposition 

i-g(c) = ^p(7) + i-c. (4) 

The cluster decomposition (4) is converged according to its definition. But 
the function Q{c) is differed from zero only at c > > and, therefore, it is not 
analytic. Some principle difficulties of the value calculation are connected 
with this surcumstance. At present time, there are not some algorithms of its 
evaluation not using the stochastic modeling. 

5. The estimate of percolation probability. For reception of the upper 
estimate of the percolation threshold, it is necessary to obtain the suitable 
estimate of the probability Q{c) from below. Such an estimate is based on the 
following statement 

L e m m a 1. It is valid the inequality 

00 

c - Q{c) < ^(1 - c)X (5) 

n=3 

where Vn = |{7 G B : I7I = n}\, n > 3. 

□ We use the elementary upper estimate P{j) < (1 — c)'"^! which follows 
from Definition 2 and the expression of ^(7). Using this estimate and (4), we 
come to the upper restriction of the sum 

00 

7eB 7eB A;=3 
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Now, we find the upper estimate of the vahie n > 3. With this aim, we 
introduce the class B„ of all such simple 0*-cycles with the length n on the 
lattice A* that each of them surrounds the point 0. Further, we introduce the 
set C„(xo) of simple ways 7 = (xq, xi..., x„) on the lattice A* having the length 
n and such that each of them possesses the property Xj 7^ j = 0, l...,n—2. 
We consider also those subsets C„_i(xo,Xi) of the last introduced set C„(o) 
which consist of ways having two leading fixed vertexes xq, Xi. Besides, each 

way of these subsets is such that each two its edges following one after another, 

00 

are necessarily a part of the cycle belonging to the family B = B^^. In 

n=3 

view of the uniformity of the hexagonal lattice (i.e. the equivalence of all its 
vertexes), the value |C„_i(xo, xi)| does not depend on the vertex xq. 

Let Q!fc = (xo,xi, ...,Xfc) is the fixed simple way such that it is the part of 
circle belonging to the class B. Further, 0(0;^) is the set of edges (xfc,X;t+i) 
such that ways ak+i = (xq, xi..., x^, x^+i) possess the same property as the 
way ak has. Denote = max{|D(Q;jfe)|}. 

L e m m a 2. It is valid the inequality 

rn<n^{n-2) max {|C„_i(xo, xi)|} . (6) 

xi:xi(^*xo 

□ Let us consider an infinite simple way a{0) from the vertex on the lattice 
A. For definiteness, this way is chosen in such a manner that each its finite part 
has the least length among all ways connecting its ends. According to Theorem 
1, each cycle 7 on A* has necessarily the vertex which is common with a{Qi). 
Among all such vertexes, we choose the vertex being the nearest one to the 
vertex where the distance is counted along the way a{Qi). We designate this 
vertex by z^. Then, each cycle 7 G B is contained in the unique class among 
not intersected classes B(/), / G N. Each class B(/) of the last family contains 
only those cycles 7 of B which have the correspondent vertex at the distance 
/ along the way q;(0), dist (z^, 0) = /. On the basis of classes B(/), / G N, we 
construct some classes B„(/) = B(/) fi B„ = {7 G B : I7I = n, dist(z^, 0) = /}, 
n = 3, 4, / G N. It is obvious that the inequality I < n — 2 takes place and. 
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n-2 

therefore, it is valid B„ = B„(/), 

1=1 

n-2 

rn = E|Bn(OI- (7) 

1=1 

Let us fix the clockwise direction on cycles 7 = (z^, xi...., x„_i, z^) of B„(/). 
We introduce into consideration disjoined sets B„(/, a) = {7 G B„(/) : xi— z^ = 
a} of cycles containing in B„(/) such that each of them has the fixed leading 
vertex that follows after z^ according to the direction on 7. This vertex is 
xi = z^ + a. The vector a is one of vectors on Fig. 2 provided that z^ is 
translated to 0. Then, the disjoint decomposition B„(/) = |^ B„(Z, y — z^) 

y:z^0*y 

takes place and, therefore, 

|B,(0|= Yl |Bn(/,y-z^)|. (8) 

It is obvious that |B„(/,a)| < |C„_i(z^,z^ + a)|. Namely, any fixed circle of 
B„(/, a) turns to the way of C„_i(z^, z^ + a) after the removal of the last circle 
edge according to the order. Thus, we obtain the injection since no more than 
one cycle of B„(/,a) corresponds to each way of C„_i(z^, z^ + a). Thus, on 
the basis of the above estimate and (8), it follows the inequality 

|B„(/)|<n,. max {|Cn-i(xo, xi)|} . 

xi:xi(j!)*xo 

Since the value |C„_i(xo, xi)| does not depend on / = dist(0, z^), we obtain 
(6) applying the received inequality to do the upper estimate of the right-hand 
side of the decomposition (7). ■ 

We introduce the construction for the estimating of the value | C„_i (xq, xi) | . 
We characterize uniquely each way of C„(xo), n > 2 by the sequence (ai, aj^_i) 
of shift vectors where ai = xi — xq and each its vector a^, i = 2..., n represents 
the vector x^ — x^ 1 turned on the angle being least among two ones formed 
by the vector Xj_i — Xi_2 and the basis vector 62 on the lattice plane. Now, we 
introduce the set G„(ai) of all sequences (ai, ...,a„) with the fixed vector ai 
and such that each pair (a^, a^+i), z = 1, 2., n — 1 is admissible in it. Here, the 
pair (a, a') is called the admissible one if there exists the circle 7 in the family 
B such that the corresponding sequence contains the pair (a, a') of vectors 
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following one after another in it. The introduced set is equivalent to the set 
C„(xo,xi). Hence, |C„(xo,Xi)| = g{3.i;n) = |G^(ai)|. 

Let us decompose the set G„(ai) by the following way G„(ai) = Gri(ai; a„) 

on not intersecting sets G,^(ai;a^) of sequences (ai, a2..., a„) where the last 
shift vector is fixed side by side the first one. Then, 5'(ai; n) = |G„(ai, a„)|, 

Introducing the numbering for possible shift vectors which is presented on 
Fig. 2, we may consider that the value (|G„(ai,a„)| = gj{ai,n);j = 1 12) is 
12-dimensional vector for each n G N. Here, the component number j = 1 12 
is defined by the number of the shift vector a„ in the accepted numbering. 

12 

Thus, g{ai;n) = ^gj{ai;n). 

Now, let us define the matrix S with the dimension equal to the number of 
nearest neighbors of the vertex on the conjugate lattice A*. For the hexagonal 
lattice, this number is equal to 12. Matrix elements Sij, i^j = 1 12 may have 
values or 1 according to the following rule. We put Sij = 1 if there is such a 
pair (afc, afc+i) of shift vectors at /c = 1, n — 1 in one of a sequence in the set 
G„(ai,a„) for anything value n = 3,4... and vectors ak,ak+i have numbers i 
and j according to the numbering accepted on Fig. 2. On the contrary, Sij = 
if there is not the circle with the external border which contains the pointed 
out joining of edges. Further, we will name the matrix S as the matrix of way 
connection. Since Sij > 0; = 1 12, then its maximal modulo eigenvalue 
Ao is positive according to the Frobenius theorem (see, [7]) for matrixes with 
nonnegative elements. 

L e m m a 3. Let the matrix S has the unique eigenvalue Aq > with 
maximal modulus. Then, the asymptotic formula 

^,(al,n) = C,,A^l(l + o(l)), n^oo (9) 

takes place where the number j corresponds to the shift vector ai and the 
nonzero matrix C has non-negative matrix elements Cij. 

□ According to the definition of the vector {gi{ai]n)] j = 1 -i- 12), the 
recurrent relation 

12 

5'j(ai;n) = ^gk{ai]n- l)Ski 
takes place for any n = 2, 3... and for the vector ai with number j. Besides, it 
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takes place ^^(ai; 1) = Sji. Then, using the induction on n G N, we conclude 
that 

p,(ai;n) = (5"-%. 

Since the matrix S has the unique eigenvalue Aq with the maximal absolute 
value, for the expression in the right-hand side of last equality, the asymptotic 
formula (9) is valid at n — ^ oo. Besides, Sij > and, hence, the nonnegativity 
of matrix elements Cij follows from the Frobenius theorem. ■ 

Now, we may find the below estimate of the probability Q{c). 

Theorems. Let the maximal eigenvalue Aq of the matrix S of way 
connections be unique. Then, it is valid the following below estimate of the 
probability Q (c) 

oo 

Qic) > c - n*(l - c)2 J2in - 2)[(1 - cjAo]""' . (10) 

1=3 

□ On the basis of Lemmas 1 and 2 and the definition of the functions 
p^(ai; n), i = 1 12, we have 

oo oo 

c - Qic) < ^(1 - c)V, < n, ^(n - 2)(1 - c)" m^ax{|C,_i(xo, xi)|} = 

n=3 n=3 ^ 

oo oo 12 

= J2in - 2)(1 - c)"^,_i(ai) = ^(n - 2)(1 - c)" ^^,(ai;n - 1) . 

n=3 n=3 i=l 

Applying the asymptotic formula (9), we obtain 

00 

c - Qic) < n,C(l - c)2 ^(n - 2)[(1 - c)Ao]"-' 

n=3 
12 



where the positive constant C > max Cij is chosen by such a way that the 
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inequality giiai^n) < CAq" , n G N takes place. ■ 

Corollary. For the percolation threshold c^ of the Bernoulli random 
field {c(x); X £ V} on the hexagonal lattice {V, the inequality c^ < 1 — Aq ^ 
is valid. 

□ The series in the right-hand side of the inequality (10) converges at 
(1 — c)Ao < 1. Then, it takes place at c > 1 — Aq ^. The convergence of 
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this series, applying the reasoning based on the Borel-Cantelli lemma (see, for 
example, pj) leads to the percolation probability Q{c) being distinct from zero. 
It takes place at the above pointed out restriction on the parameter c. Hence, 
< 1 - Aq ^ ■ 

6. The upper estimate of the percolation threshold. It is follows 
from considerations of the previous section that the problem of the upper 
estimation of the percolation threshold on the hexagonal lattice is reduced to 
calculation of the eigenvalue Aq. For its calculation, first of all, it is necessary 
to find the matrix S. For the evaluation of matrix elements Sij^ it is necessary 
either to find the cluster which has the external border containing the pair 
of edges (a, a') where vectors a, a' have numbers i and j correspondingly or 
to prove that there is not the cluster with such a pair. Following statements 
describe the general structure of the matrix S and point out its zero elements. 



L e m m a 4. The matrix Sij is symmetric and Sa = 0. 

□ Since the passage of the cycle being the external border is possible in 
both directions, the existence or the absence of pair edges (a, a') following 
one after another in it which have numbers i and j correspondingly leads to 
existence or absence of the cycle which has the pair of edges (a', a). It means 

The equality Sa = 1 corresponds to the fact that there is the pair of edges 
in the external border which are described by the pair (a, a) of shift vectors 
in the sequence (ai..., a„) and the vector a has number i. It means, according 
to the description of the cycle by means of the specified sequence that these 
edges look as (x, y) and (y,x). But it is impossible. Therefore, Sa = for all 
values I = 1^12. ■ 

L e m m a 5. The Matrix Sij has the property Sij = *S'i^i4_j, S^j = S^^iQ-j, 
Sgj = Sg^is^j, j = 2 12 where the subtraction operation in bottom indexes 
is understood modulo 12. 

□ The first relation follows from the symmetry of Fig. 2 describing the join- 
ing of edges in the vertex relative to the reflection from the straight line 
defined by vertexes 1, 0, 7. The second relation follows from the symmetry 
of this figure to the reflection from the straight line passing through vertexes 
5, 0, 11 and the third relation is connected with the reflection from the straight 
line passing through vertexes 3,0,9. ■ 

L e m m a 6. Matrix elements Sij have the property 5'i-|_4j+4 = Sij where 
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sums 2 + 4 and j + 4: are understood modulo 12. 

□ This property follows from the symmetry of Fig. 2 describing all possible 
connections of edges in the vertex relative to rotations on the angle 27r/3. ■ 




Corollary. The matrix S has the foUowing block structure made of 
two 4 X 4-matrices F, G 

f F G G+\ 
S= G+ F G . (11) 
\G G+ F J 

□ Using 4 X 4-matrices F, G, H with matrix elements Fij = Sij, Gij = Sij+4, 
Hij = Sij+8, hj = 1-^4, it follows from Lemma 6 that the matrix S is 
represented as 

/ F G H 
S= H F 
\G H 

From Lemma 4, it follows that H = G'^, as Gij = = Si+gj si i,j = 1 -i- 4 
and, simultaneously, Si+s,j = Sj^i+s = H'^- B 
L e m m a 7. Si^^ = 0. 

□ Each vertex x of the external border dW should have the vertex y adja- 
cent with it and belonging to W. Besides, there exists the vertex z adjacent 
with X, but not belonging W U dW. If we admit that 5i,5 = 1, then vertexes 
1 and 5 (see. Fig. 3 where joinings of edges in the vertex are shown by the 
dotted line) belong to dW and, hence, they do not belong to W. Since the 
vertex has three vertexes 1,5,9 adjacent with it, then the vertex 9 should be 
simultaneously the vertex of W and it does not belong to U dW. But, it is 
impossible. ■ 

L e m m a 8. Sij = 0, i, j = 1 4. 

□ It is sufficient to prove that Sij = 0, j = 2, 3, 4 and S2,j = 0, j = 3, 4. We 
consider the first group of matrix elements. Vertexes and 1 belong to dW. 
There are vertexes x, y G such that x01, y^O and also there are vertexes 
x.' ,-y' ^ W U dW such that x'0 1, y'cj) 0. Further, there are ways from vertexes 
x' and y' to infinity not crossing W. Vertexes 2 and 12 may play the role of 
vertexes x, x' for the vertex 1 and only vertexes 5 and 9 may be the same for 
the vertex 0. Then, there are four choice variants of vertex pairs x, x' and y, y' 
among presented possibilities. It is sufficient to analyze only two cases which 
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Figure 3: Possible joinings of edges in the vertex 0. 

do not reduce to each other by mirror reflection relative to the straight hne 
passing through vertexes 1,0,7. We choose the foflowing possibihties. Vertexes 
2 and 12 are the vertexes x, x' in both cases. And the vertexes 9 and 5 are 
the vertexes y, y' in the flrst case and, on the contrary, they are interchanged 
their position in the second. 

In the flrst case, there are not inflnite ways from vertexes 1 and which are 
crossed with W and passed through vertexes 12 and 5 correspondingly. We 
consider the line in the plane where the lattice is placed. This line consists of 
the flrst inflnite way of lattice edges that comes from inflnity to the vertex 12. 
Further, it consists of the sequence of edges (12, 1), (1, 0), (0, 5) and it is ended 
by the second inflnite way from the vertex 5 to inflnity. The constructed line 
divides the plane into two parts so that vertexes x and y are in different ones. 
Then, vertexes x and y should not belong to the same cluster W since any way 
from the vertex x to the vertex y on lattice edges crosses necessarily this line. 
Besides, it may be done only in the lattice vertex. Hence, the edge connection 
(1,0) with one of edges (0,j), j = 2,3,4 is impossible in the variant under 
consideration. 

In the second case, it is considered the analogous line on the plane which 
consists of the way coming from inflnity to the vertex 1 without crossing of 
the cluster W. Further, it consists of the edge (1,0), the diagonal of the given 
lattice hexagon which is led from the vertex to the vertex j (it is presented 
on Fig.3 where this side has the border vertexes 0, 1, 2, 3, 4, 5) and the way 
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from the vertex j to infinity on lattice edges without the crossing of W. The 
fine constructed divides the plane into two parts so that vertexes x and y are 
in diff^erent ones. Therefore, they should not belong to the same cluster W 
since any way on lattice edges from the vertex y to the vertex x necessarily 
crosses the line in the lattice vertex. 

Contradictions obtained in both cases show that Si^j = at j = 2, 3, 4. 

Let us consider the second group of matrix elements 5*23, S24. We construct 
the following line on the lattice plane. Firstly, it consists of the way which 
comes to the vertex 2 from infinity without crossing of the cluster W. It 
should exist since the vertex 2 belongs to dW. Further, the line consists of 
the consecutive passage of two diagonals on the hexagon presented on Fig. 3. 
Diagonals are defined by vertex pairs (2, 0) and (0, j) where j = 3, 4. The line 
is ended by the way from the vertex j to infinity. It divides the plane into two 
parts and the cluster W should be settled down completely in one of them. 

It is obvious that 5*23 = since the cluster W may be not settled down on 
the right-hand side of the constructed line according to the direction accepted 
on it. Otherwise, the vertex has no vertexes belonging to W and being 
adjacent with it. On the other hand, the cluster W may not be on the left- 
hand side of this line since there is the way leaving from the vertex to infinity 
without crossing of the cluster W. It should exist according to the definition 
of the external border vertex. This way should be also in the left part of the 
plane. Then, the way divides the left-hand side into two parts again. The 
cluster W should be settled down completely in one of them and, therefore, 
one of vertexes 2 or j does not belong to the external border since there is not 
anything vertex of W adjacent to it. The obtained contradiction proves the 
equality ^24 = O.B 

Theorem4. The matrix S is presented by the formula 




S = \ G+ G , G = 



(0 


1 


1 


1\ 





1 


1 


1 





1 


1 


1 




1 


1 





(12) 



where is the zero matrix. 

□ It follows from Lemmas 5 and 8 that Ss,a = Ss,2 = 82,3 = 0. Then the 
matrix F in the formula (11) is zero owing to Lemma 4. On the basis of the 
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property Sij = 5'i,i4_j and Lemmas 7 and 8, we find 

'S'1,12 = 'S'1,2 = 0; ^i^ii = 51,3 = 0; .Si^io = »S'i,4 = 0; ^i^g = 51,5 = . 

Then, the first Une in the matrix G'^ in the formula (11) and, hence, the first 
column in the matrix G are zero. Other elements of the matrix G are equal 
to unity due to the above mentioned criterion at the evaluation of matrix 
elements. So, it is had Sij = 1 ai i = 1 4,j = 6, 7, 8 (see Figs. 4-7). ■ 




Figure 4: Clusters containing the edge 1-0. 




Figure 5: Clusters containing the edge 2-0. 



Notice. It follows from the matrix S expression that = 7. 
On the basis of Theorem II in Appendix, having put n = 3, m = 4, ^4 = <S, 

we conclude that the maximal eigenvalue Aq of 
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Figure 6: Clusters containing the edge 3-0. 




Figure 7: Clusters containing the edge 4-0- 



the matrix S coincides with the maximal eigenvalue of the matrix 



5 = + G' + G'+ = 



/O 

1 
1 



1 

2 
2 
2 



1 
2 
2 
2 



1\ 

2 
2 

2/ 



The rank of this matrix is obviously equal 2. Then, it has two zero eigenvalues 
and, therefore, its characteristic equation looks as follows 

det(B - Al) = (A^ - A^i + ^2)A2 = 

where = Sp5, .^2 = [Spfi^ - [S^Bf]/2. Then, = Sp5 = 6 and 6 = 3 
since diagonal elements of the matrix 5^ are equal (B^)^j = 3, 13, 13, 13. Hence, 
the number Aq is the greatest root of the quadratic equation A^ — 6A — 3 = 
0. Whence, we obtain Aq = 3 + 2\/3. Thus, we have proved the following 
statement. 
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Theorems. The maximal eigenvalue Xq of the matrix S is equal 
3 + 2^3^6,46. 

From Theorem 5 and Corollary of Theorem 3, it follows directly 
The basic statement. The percolation threshold of the Bernoulli field 
on the hexagonal lattice does not surpass the number 2/3(3 — v^). 
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Appendix 

L e m m a I. Let B he the m x m-matrix with nonnegative elements such 
that there is the number /j, > and the vector b G M™, b = (61, ...,bm) with 
nonnegative components such that it takes place the inequality 

m 

1=1 

Then the maximal eigenvalue X{B) of the matrix B on the absolute value is 
positive and it satisfies the inequality X{B) > ji. 

□ We apply N times the inequality (Al). Using the induction on A/", it is 
obtained the inequality 

m 

1=1 

At first, we prove the inequality (Al) in the case when the matrix B has only 
one-dimensional eigenspaces. Let c^'"-' = {c^i^^]k = 1 m) be eigenvectors of 
the matrix B with correspondent eigenvalues A^., r = 1 -i- m. We consider that 
all Xr have been ordered according to the decrease of their absolute values (at 
the coincidence of absolute values, they are ordered according to their phases). 
We write down the decomposition of the vector 

m 
r=l 

Then, 

m 
r=l 

and, therefore, the inequality (A2) is represented in the form 

TO 

M%<E«-^M''- (^3) 

r=l 

There is the number k < m such that 6^ > 0. Let us consider the inequality 
(A3) for this number. Among all numbers r = 1 m, there is the minimal 
number tq such that ^roC^^°'' 7^ 0- Then, the inequality (A3) is rewritten in the 
form 

TO 

r=ro 
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or, calculating the N-th degree root for both positive parts of the inequality, 
we obtain 



k 



Lr=ro 



1/N 



= lA, 



k 



r=ro 



1/N 



The last equality is valid due to the positivity of the right-hand side of (A3). 
Further, we go to the limit N ^ oo. In this case, lim bV^ = 1. Due to the 

boundedness of the summing expression in the right-hand side of the inequality, 
its upper limit does not surpass the unity. Then, /i < \Xro\- According to the 
definition and due to the Frobenius theorem, |Arp| < X{B). 

One may change any matrix B by such a way that all its eigenspaces turns 
into one-dimensional ones. It is done by means of addition to the matrix B of 
the suitable matrix with nonnegative elements. It may be done as much near 

2 

to zero matrix as one wants. This nearness is understood according to the MJ^ 

topology. The possibility of this follows from the fact that the multiplicity 

condition of eigenvalues satisfying the characteristic equation det{B — XI) = 

d 

is expressed by the supplement equation — det{B — XI) = 0. The last equation 



dX 



cuts out a differential manifold with the codimensionality one in the space 
of admissible matrices B. At such an addition of small matrix, the matrix B 
transforms to such a matrix B^ that 5g ^ 5 at e — ^ +0. 

The inequality (Al) is proved in general case by the following way. Let b 
be the eigenvector of the matrix B corresponding to X{B). Since 



m 



XiB)bk = J2Bkibi<J2(Be)kibi 



1=1 



1=1 



and the matrix B^ has only one-dimensional eigenspaces, the inequality (Al) 
takes place where it is needed to change both X{B) on X{B^) and /i on X{B). 
Consequently, X{B) < X{B^). After that, we come to the limit e — ^ +0. ■ 

Theoreml. Let A be the n x n-matrix with nonnegative elements. It 
consists of w? blocks being rectangular matrixes B^^'^^ which have correspond- 
ingly pk lines and si columns k,l = 1 m, m < n, pi -\- p2 + ... + Pm 



n, 



20 



Si + S2 + •• 



= n, 



A = 



I 5(1,1) ^(1,2) 
5(2,1) ^(2,2) 



5(1,H \ 
^(2,m) 



\ 



(^4) 



Matrix elements of the square m x m-matrix B are defined by tlie formula 



j=i 



Pk 



/c, / = 1 m 



Tlien, positive eigenvalues X{A) and X{B) corresponding to matrixes A and B 
which are maximal on their absolute values satisfy the inequality X{B) > X{A) . 

□ Matrixes A and B possess nonnegative elements. Due to the Frobenius 
theorem jTj, eigenvalues of these matrixes are positive if they are maximal on 
their absolute values. Besides, due to this theorem, there are such eigenvectors 
of matrixes A and B corresponding to eigenvalues pointed out that they have 
nonnegative components. We designate them as follows a = (ai, a„), aj > 0, 
j = 1 n; b = (61, &to), 6i > 0, i = 1 m. Then, 



j=i 



m 



X{A)a,, J2Bkibi = X{B)bk 



1=1 



We introduce the vector c G M"^, c = (ci,...,Cto) with components Ck = 
maxja^; si + ... + Sfc_i < ^ < si + ... + s/j}. Dividing summations on j = in 
both sides of first equality, we introduce repeated summations. First of them 
is done on groups containing si, numbers, the second is done within each 
of these groups. In a result, we obtain 



j = l j = Si + ... + S;_i + l / = 1 j = Si 



S1 + ... + S1 

E 



Using the change i =^ pi + ...+pk-i + i, i = l-^Pk^ we rewrite the last inequality 
in the form 

m si 

X{A)a,+p^+...+p,_^ < ^Q^bI 

1=1 j=i 
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Calculating the maximum on i in both sides of the inequality, we obtain 



A(A)cfc < V Bk, 



1=1 



ici 



(A5) 



The last inequality coincides with (Al), if we put /i = X{A) and change b on 
c. Then, the theorem statement follows from (Al). ■ 

T h e o r e m II. Let {A^^\ ...^A^^'') be the ordered collection of m x m- 
matrixes with nonnegative elements. Further, let nm x nm-matrix A be made 
up of m X m-matrixes B^'^'^^ according to the following formula 



A = 



5(2,1) 



5(1,2) 
5(2,2) 



5(l,n) \ 
5(2,n) 



V 



^(n,l) 5('^,2) 



(AQ) 



I 



where each of ordered collections {B^'^'^\ 5*^*'"^), i = 1, ...,n is obtained by 



the permutation P G of the collection {A^^\ A^"^^) and P„ is the permu- 
tation group of the n-th order. Then, the maximal eigenvalue of the matrix A 
coincides with the maximal eigenvalue of the n x n-matrix B = A^-*^^ + . . . + A^'^'' . 



□ The matrix A has nonnegative elements. Therefore, according to the 
Frobenius theorem, its eigenvalue X{A) being maximal on the absolute value is 
positive. Let a G M"™, a = (ai, ...,anm) be the eigenvector which corresponds 
to this eigenvalue where > 0, z = 1 ^ nm, Aa = X{A)a. It means that 

nm 

Aijaj = X{A)ai , i = 1 nm . 

i=i 

We define the vector c G W^, c = (ci, c^), cj > 0, j = 1 m where 

Ck = max{afc+sm; s = 0, ...,n-l}, k = l^m. 
Then, changing i = k -\- pm, j = 1 + sm, s,p = 0, 1, n — 1, we have 

1 m m / n—l \ 

X{A)ak+pm — ^ ^ ^ ^ ^A;+pm,Z+smQ'/+5m ^ ^ ^ { ^ ^ ^fc+pm,Z+sm 1 Q • 
s=0 Z=l Z=l \s=0 / 
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Calculating the maximum on p = 0, 1, . . . , n — 1 of both sides of the inequality, 
we obtain 

A(A)cfc < Bkici 
1=1 

where 

Bki = max |y^^fc+w+sm;P = 0, 1, ...,n- l| {A7) 

are matrix elements B = A^^^ + ... + A^'^\ This formula follows from the 
equality 

n—l n—1 
s=0 s=0 

where the permutation P G P„ is defined by the number p being the line 
number in the block matrix ( A6) . Having changed the summation variable as 
follows P(s + 1) ^ s, we find that the last expression is equal to 



n—l n 



s=0 s=l 



and the sum in (A4) does not depend on p. 

Applying the statement of Theorem I to the matrix A, in a result, we obtain 
the inequality X{B) > X{A). On the other side, we take the eigenvector b of 
the matrix B with nonnegative components that corresponds to the eigenvalue 
X{B). Further, we define the vector c' = (b, b, . . . , b). For this vector, we have 

n 

Ac' = (Bb, Bb, . . . , Bb) = X{B)c' . 

Then, X{B) is the eigenvalue of the matrix A with the eigenvector c'. According 
to the definition of eigenvalue X{A), we have X{A) > X{B). From two obtained 
inequalities, it follows that X{A) = X{B). ■ 
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